Experimental quantification and modelling of attrition of infant formulae during pneumatic conveying

Infant formula is often produced as an agglomerated powder using a spray drying process. Pneumatic conveying is commonly used for transporting this product within a manufacturing plant. The transient mechanical loads imposed by this process cause some of the agglomerates to disintegrate, which has implications for key quality characteristics of the formula including bulk density and wettability. This thesis used both experimental and modelling approaches to investigate this breakage during conveying. One set of conveying trials had the objective of establishing relationships between the geometry and operating conditions of the conveying system and the resulting changes in bulk properties of the infant formula upon conveying. A modular stainless steel pneumatic conveying rig was constructed for these trials. The mode of conveying and air velocity had a statistically-significant effect on bulk density at a 95% level, while mode of conveying was the only factor which significantly influenced D[4,3] or wettability. A separate set of conveying experiments investigated the effect of infant formula composition, rather than the pneumatic conveying parameters, and also assessed the relationships between the mechanical responses of individual agglomerates of four infant formulae and their compositions. The bulk densities before conveying, and the forces and strains at failure of individual agglomerates, were related to the protein content. The force at failure and stiffness of individual agglomerates were strongly correlated, and generally increased with increasing protein to fat ratio while the strain at failure decreased. Two models of breakage were developed at different scales; the first was a detailed discrete element model of a single agglomerate. This was calibrated using a novel approach based on Taguchi methods which was shown to have considerable advantages over basic parameter studies which are widely used. The data obtained using this model compared well to experimental results for quasi-static uniaxial compression of individual agglomerates. The model also gave adequate results for dynamic loading simulations. A probabilistic model of pneumatic conveying was also developed; this was suitable for predicting breakage in large populations of agglomerates and was highly versatile: parts of the model could easily be substituted by the researcher according to their specific requirements.

Applications of finite element computer modelling and thermal infrared remote sensing to the study of geothermal anomalies

Buried heat sources can be investigated by examining thermal infrared images and comparing these with the results of theoretical models which predict the thermal anomaly a given heat source may generate. Key factors influencing surface temperature include the geometry and temperature of the heat source, the surface meteorological environment, and the thermal conductivity and anisotropy of the rock. In general, a geothermal heat flux of greater than 2% of solar insolation is required to produce a detectable thermal anomaly in a thermal infrared image. A heat source of, for example, 2-300K greater than the average surface temperature must be a t depth shallower than 50m for the detection of the anomaly in a thermal infrared image, for typical terrestrial conditions. Atmospheric factors are of critical importance. While the mean atmospheric temperature has little significance, the convection is a dominant factor, and can act to swamp the thermal signature entirely. Given a steady state heat source that produces a detectable thermal anomaly, it is possible to loosely constrain the physical properties of the heat source and surrounding rock, using the surface thermal anomaly as a basis. The success of this technique is highly dependent on the degree to which the physical properties of the host rock are known. Important parameters include the surface thermal properties and thermal conductivity of the rock. Modelling of transient thermal situations was carried out, to assess the effect of time dependant thermal fluxes. One-dimensional finite element models can be readily and accurately applied to the investigation of diurnal heat flow, as with thermal inertia models. Diurnal thermal models of environments on Earth, the Moon and Mars were carried out using finite elements and found to be consistent with published measurements. The heat flow from an injection of hot lava into a near surface lava tube was considered. While this approach was useful for study, and long term monitoring in inhospitable areas, it was found to have little hazard warning utility, as the time taken for the thermal energy to propagate to the surface in dry rock (several months) in very long. The resolution of the thermal infrared imaging system is an important factor. Presently available satellite based systems such as Landsat (resolution of 120m) are inadequate for detailed study of geothermal anomalies. Airborne systems, such as TIMS (variable resolution of 3-6m) are much more useful for discriminating small buried heat sources. Planned improvements in the resolution of satellite based systems will broaden the potential for application of the techniques developed in this thesis. It is important to note, however, that adequate spatial resolution is a necessary but not sufficient condition for successful application of these techniques.

Studies relating to the population dynamics of the European badger (Meles meles) in Ireland

European badgers (Meles meles) are an important part of the Irish ecosystem; they are a component of Ireland’s native fauna and are afforded protection by national and international laws. The species is also a reservoir host for bovine tuberculosis (bTB) and implicated in the epidemiology of bTB in cattle. Due to this latter point, badgers have been culled in the Republic of Ireland (ROI) in areas where persistent cattle bTB outbreaks exist. The population dynamics of badgers are therefore of great pure and applied interest. The studies within this thesis used large datasets and a number of analytical approaches to uncover essential elements of badger populations in the ROI. Furthermore, a review and meta-analysis of all available data on Irish badgers was completed to give a framework from which key knowledge gaps and future directions could be identified (Chapter 1). One main finding suggested that badger densities are significantly reduced in areas of repeated culling, as revealed through declining trends in signs of activity (Chapter 2) and capture numbers (Chapter 2 and Chapter 3). Despite this, the trappability of badgers was shown to be lower than previously thought. This indicates that management programmes would require repeated long-term efforts to be effective (Chapter 4). Mark-recapture modelling of a population (sample area: 755km2) suggested that mean badger density was typical of continental European populations, but substantially lower than British populations (Chapter 4). Badger movement patterns indicated that most of the population exhibited site fidelity. Long-distance movements were also recorded, the longest of which (20.1km) was the greatest displacement of an Irish badger currently known (Chapter 5). The studies presented in this thesis allows for the development of more robust models of the badger population at national scales (see Future Directions). Through the use of large-scale datasets future models will facilitate informed sustainable planning for disease control.

Grid integration of wave energy & generic modelling of ocean devices for power system studies

The work presented in this thesis covers four major topics of research related to the grid integration of wave energy. More specifically, the grid impact of a wave farm on the power quality of its local network is investigated. Two estimation methods were developed regarding the flicker level Pst generated by a wave farm in relation to its rated power as well as in relation to the impedance angle ψk of the node in the grid to which it is connected. The electrical design of a typical wave farm design is also studied in terms of minimum rating for three types of costly pieces of equipment, namely the VAr compensator, the submarine cables and the overhead line. The power losses dissipated within the farm's electrical network are also evaluated. The feasibility of transforming a test site into a commercial site of greater rated power is investigated from the perspective of power quality and of cables and overhead line thermal loading. Finally, the generic modelling of ocean devices, referring here to both wave and tidal current devices, is investigated.

Preparation and characterisation of ceria particles

Cerium dioxide (ceria) nanoparticles have been the subject of intense academic and industrial interest. Ceria has a host of applications but academic interest largely stems from their use in the modern automotive catalyst but it is also of interest because of many other application areas notably as the abrasive in chemical-mechanical planarisation of silicon substrates. Recently, ceria has been the focus of research investigating health effects of nanoparticles. Importantly, the role of non-stoichiometry in ceria nanoparticles is implicated in their biochemistry. Ceria has well understood non-stoichiometry based around the ease of formation of anion vacancies and these can form ordered superstructures based around the fluorite lattice structure exhibited by ceria. The anion vacancies are associated with localised or small polaron states formed by the electrons that remain after oxygen desorption. In simple terms these electrons combine with Ce4+ states to form Ce3+ states whose larger ionic radii is associated with a lattice expansion compared to stoichiometric CeO2. This is a very simplistic explanation and greater defect chemistry complexity is suggested by more recent work. Various authors have shown that vacancies are mobile and may result in vacancy clustering. Ceria nanoparticles are of particular interest because of the high activity and surface area of small particulates. The sensitivity of the cerium electronic band structure to environment would suggest that changes in the properties of ceria particles at nanoscale dimensions might be expected. Notably many authors report a lattice expansion with reducing particle size (largely confined to sub-10 nm particles). Most authors assign increased lattice dimensions to the presence of a surface stable Ce2O3 type layer at low nanoparticle dimensions. However, our understanding of oxide nanoparticles is limited and their full and quantitative characterisation offers serious challenges. In a series of chemical preparations by ourselves we see little evidence of a consistent model emerging to explain lattice parameter changes with nanoparticle size. Based on these results and a review of the literature it is worthwhile asking if a model of surface enhanced defect concentration is consistent with known cerium/cerium oxide chemistries, whether this is applicable to a range of different synthesis methods and if a more consistent description is possible. In Chapter one the science of cerium oxide is outlined including the crystal structure, defect chemistry and different oxidation states available. The uses and applications of cerium oxide are also discussed as well as modelling of the lattice parameter and the doping of the ceria lattice. Chapter two describes both the synthesis techniques and the analytical methods employed to execute this research. Chapter three focuses on high surface area ceria nano-particles and how these have been prepared using a citrate sol-gel precipitation method. Changes to the particle size have been made by calcining the ceria powders at different temperatures. X-ray diffraction methods were used to determine their lattice parameters. The particles sizes were also assessed using transmission electron microscopy (TEM), scanning electron microscopy (SEM), and BET, and, the lattice parameter was found to decrease with decreasing particle size. The results are discussed in light of the role played by surface tension effects. Chapter four describes the morphological and structural characterization of crystalline CeO2 nanoparticles prepared by forward and reverse precipitation techniques and compares these by powder x-ray diffraction (PXRD), nitrogen adsorption (BET) and high resolution transmission electron microscopy (HRTEM) analysis. The two routes give quite different materials although in both cases the products are essentially highly crystalline, dense particulates. It was found that the reverse precipitation technique gave the smallest crystallites with the narrowest size dispersion. This route also gave as-synthesised materials with higher surface areas. HRTEM confirmed the observations made from PXRD data and showed that the two methods resulted in quite different morphologies and surface chemistries. The forward route gives products with significantly greater densities of Ce3+ species compared to the reverse route. Data are explained using known precipitation chemistry and kinetic effects. Chapter five centres on the addition of terbia to ceria and has been investigated using XRD, XRF, XPS and TEM. Good solid solutions were formed across the entire composition range and there was no evidence for the formation of mixed phases or surface segregation over either the composition or temperature range investigated. Both Tb3+ and Tb4+ ions exist within the solution and the ratios of these cations are consistent with the addition of Tb8O15 to the fluorite ceria structure across a wide range of compositions. Local regions of anion vacancy ordering may be visible for small crystallites. There is no evidence of significant Ce3+ ion concentrations formed at the surface or in the bulk by the addition of terbia. The lattice parameter of these materials was seen to decrease with decreasing crystallite size. This is consistent with increased surface tension effects at small dimension. Chapter six reviews size related lattice parameter changes and surface defects in ceria nanocrystals. Ceria (CeO2) has many important applications, notably in catalysis. Many of its uses rely on generating nanodimensioned particles. Ceria has important redox chemistry where Ce4+ cations can be reversibly reduced to Ce3+ cations and associated anion vacancies. The significantly larger size of Ce3+ (compared with Ce4+) has been shown to result in lattice expansion. Many authors have observed lattice expansion in nanodimensioned crystals (nanocrystals), and these have been attributed to the presence of stabilized Ce3+ -anion vacancy combinations in these systems. Experimental results presented here show (i) that significant, but complex changes in the lattice parameter with size can occur in 2-500 nm crystallites, (ii) that there is a definitive relationship between defect chemistry and the lattice parameter in ceria nanocrystals, and (iii) that the stabilizing mechanism for the Ce3+ -anion vacancy defects at the surface of ceria nanocrystals is determined by the size, the surface status, and the analysis conditions. In this work, both lattice expansion and a more unusual lattice contraction in ultrafine nanocrystals are observed. The lattice deformations seen can be defined as a function of both the anion vacancy (hydroxyl) concentration in the nanocrystal and the intensity of the additional pressure imposed by the surface tension on the crystal. The expansion of lattice parameters in ceria nanocrystals is attributed to a number of factors, most notably, the presence of any hydroxyl moieties in the materials. Thus, a very careful understanding of the synthesis combined with characterization is required to understand the surface chemistry of ceria nanocrystals.

Bottom-up modelling of energy demand and technical energy savings potential in the Irish residential sector

The International Energy Agency has repeatedly identified increased end-use energy efficiency as the quickest, least costly method of green house gas mitigation, most recently in the 2012 World Energy Outlook, and urges all governing bodies to increase efforts to promote energy efficiency policies and technologies. The residential sector is recognised as a major potential source of cost effective energy efficiency gains. Within the EU this relative importance can be seen from a review of the National Energy Efficiency Action Plans (NEEAP) submitted by member states, which in all cases place a large emphasis on the residential sector. This is particularly true for Ireland whose residential sector has historically had higher energy consumption and CO2 emissions than the EU average and whose first NEEAP targeted 44% of the energy savings to be achieved in 2020 from this sector. This thesis develops a bottom-up engineering archetype modelling approach to analyse the Irish residential sector and to estimate the technical energy savings potential of a number of policy measures. First, a model of space and water heating energy demand for new dwellings is built and used to estimate the technical energy savings potential due to the introduction of the 2008 and 2010 changes to part L of the building regulations governing energy efficiency in new dwellings. Next, the author makes use of a valuable new dataset of Building Energy Rating (BER) survey results to first characterise the highly heterogeneous stock of existing dwellings, and then to estimate the technical energy savings potential of an ambitious national retrofit programme targeting up to 1 million residential dwellings. This thesis also presents work carried out by the author as part of a collaboration to produce a bottom-up, multi-sector LEAP model for Ireland. Overall this work highlights the challenges faced in successfully implementing both sets of policy measures. It points to the wide potential range of final savings possible from particular policy measures and the resulting high degree of uncertainty as to whether particular targets will be met and identifies the key factors on which the success of these policies will depend. It makes recommendations on further modelling work and on the improvements necessary in the data available to researchers and policy makers alike in order to develop increasingly sophisticated residential energy demand models and better inform policy.

Layered graphitic carbon host formation during liquid-free solid state growth of metal pyrophosphates

We report a successful ligand- and liquid-free solid state route to form metal pyrophosphates within a layered graphitic carbon matrix through a single step approach involving pyrolysis of previously synthesized organometallic derivatives of a cyclotriphosphazene. In this case, we show how single crystal Mn2P2O7 can be formed on either the micro- or the nanoscale in the complete absence of solvents or solutions by an efficient combustion process using rationally designed macromolecular trimer precursors, and present evidence and a mechanism for layered graphite host formation. Using in situ Raman spectroscopy, infrared spectroscopy, X-ray diffraction, high resolution electron microscopy, thermogravimetric and differential scanning calorimetric analysis, and near-edge X-ray absorption fine structure examination, we monitor the formation process of a layered, graphitic carbon in the matrix. The identification of thermally and electrically conductive graphitic carbon host formation is important for the further development of this general ligand-free synthetic approach for inorganic nanocrystal growth in the solid state, and can be extended to form a range of transition metals pyrophosphates. For important energy storage applications, the method gives the ability to form oxide and (pyro)phosphates within a conductive, intercalation possible, graphitic carbon as host–guest composites directly on substrates for high rate Li-ion battery and emerging alternative positive electrode materials

Modelling and predictive control techniques for building heating systems

Model predictive control (MPC) has often been referred to in literature as a potential method for more efficient control of building heating systems. Though a significant performance improvement can be achieved with an MPC strategy, the complexity introduced to the commissioning of the system is often prohibitive. Models are required which can capture the thermodynamic properties of the building with sufficient accuracy for meaningful predictions to be made. Furthermore, a large number of tuning weights may need to be determined to achieve a desired performance. For MPC to become a practicable alternative, these issues must be addressed. Acknowledging the impact of the external environment as well as the interaction of occupants on the thermal behaviour of the building, in this work, techniques have been developed for deriving building models from data in which large, unmeasured disturbances are present. A spatio-temporal filtering process was introduced to determine estimates of the disturbances from measured data, which were then incorporated with metaheuristic search techniques to derive high-order simulation models, capable of replicating the thermal dynamics of a building. While a high-order simulation model allowed for control strategies to be analysed and compared, low-order models were required for use within the MPC strategy itself. The disturbance estimation techniques were adapted for use with system-identification methods to derive such models. MPC formulations were then derived to enable a more straightforward commissioning process and implemented in a validated simulation platform. A prioritised-objective strategy was developed which allowed for the tuning parameters typically associated with an MPC cost function to be omitted from the formulation by separation of the conflicting requirements of comfort satisfaction and energy reduction within a lexicographic framework. The improved ability of the formulation to be set-up and reconfigured in faulted conditions was shown.